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The stochastic dynamics of coUoidal particles with surface activity — in the form of catalytic re- 
action or particle release — and self-phoretic effects is studied analytically. Three different time 
scales corresponding to inertial effects, solute redistribution, and rotational diffusion are identified 
and shown to lead to a plethora of different regimes involving inertial, propulsive, anomalous, and 
diffusive behaviors. For symmetric active colloids, a regime is found where the mean-squared dis- 
placement has a super-diffusive t^^^ behavior. At the longest time scales, an effective diffusion 
coefficient is found which has a non-monotonic dependence on the size of the colloid. 
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The development of biomimetic technology would sig- 
nificantly benefit from the ability to make synthetic com- 
ponents with desired motility properties. In recent years, 
there has been a range of developments along these lines, 
with functionalities that can be manipulated at length 
scales ranging from microns down to molecular scales 
These include experimental realization of actuated 
microswimmers Q and theoretical proposals of simple 
model swimmers that can tackle the problems caused 
by low Reynolds number conditions [3|. Moreover, it 
has been recently demonstrated (both experimentally 
and theoretically) that intcrfacial phoretic effects (such 
as electrophoresis, clcctroosmosis, and diffusiophoresis) 
could lead to self-propulsion of colloidal particles [1, H, 3] ■ 
It has also been shown that phoretic effects can be used 
to steer both active and passive colloidal particle 0, [1] , 
which adds to the promise of these phenomena for de- 
signing functional sclf-motilc vessels in the nanoscale. 

A fundamental property of such small objects, even 
when equipped with a self-propulsion mechanism, is that 
their motion is stochastic due to the ambient fluctua- 
tions that could be of thermal origin or otherwise. This 
means that we cannot directly control the motion of self- 
motile objects, and any design characteristic needs to be 
incorporated into statistical average outcomes. For ex- 
ample, it has been shown that a self-propelled colloidal 
particle makes a crossover between ballistic and diffusive 
behaviors over a time scale that is set by the rotational 
diffusion of the colloid, when its orientation will be ran- 
domized 0. A study of a two dimensional model of self- 
propelled objects with fluctuations both in direction and 
magnitude of the velocity has shown the possibility of 
re-entrant ballistic and diffusive behaviors [9| . Other ex- 
amples include the effect of the activity of proteins on 
the dynamics of membranes [lo| . and collective behav- 
ior of mixtures of motors and filaments (or reorganizing 
liviiig cells) and active particles such as swimming bacte- 
ria [ill . [T3 | , where a host of qualitative and quantitative 
changes have been found to occur due to nonequilibrium 
fluctuations. In light of this inherent feature, it will be 
natural to ask how many distinct regimes of motion we 



could have for active colloidal particles, what the relevant 
time scales that differentiate between these regimes are, 
and how they can be tuned so that the desired type of 
motion could be achieved by choosing the right parame- 
ters. 

Here, we aim to address some of these questions for a 
class of isolated self-motile active colloids. We study the 
velocity autocorrelation function and the mean-squared 
displacement of surface-active spherical colloidal parti- 
cles that interact with their self-generated surrounding 
clouds of solute particles via interfacial phoretic effects. 
We identify the relevant time scales in the dynamics, 
namely the hydrodynamic relaxation time that con- 
trols the crossover between inertial and viscous regimes, 
the diffusion time of the solute particles around the col- 
loid Trf, and the rotational diffusion time of the colloid r^. 
We calculate the contribution due to hydrodynamic fluc- 
tuations, as well as the self-phoretic contributions that 
depend on whether the particles are symmetric (in which 
case there is no net propulsion) or asymmetric (where the 
colloids are self-propelled). We flnd that these different 
contributions lead to a variety of different regimes, as 
summarized in Fig. [1] For symmetric surface activity, 
we flnd a regime corresponding to <C t <C where 
the active colloid demonstrates a super-diffusive behav- 
ior with a mean-squared displacement ~ t^/^ (shaded 
region in Fig. [Ij. This is similar and somewhat related 
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FIG. 1: Summary of results for the different contributions to 
the mean-squared displacement of active colloids. The total 
mean-squared displacement is obtained by summing all of the 
contributions for asymmetric colloids, and the bottom two 
rows (only) for symmetric colloids. 
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FIG. 2: (a) Schematics of an instantaneous configuration of 
an axially symmetric surface-active spherical colloid, (b) The 
model could correspond to a chemical reaction catalyzed on 
the surface with the simplifying assumption that one of the 
product particles is very similar to the substrate, or a con- 
tainer that releases particles through channels. 

to an anomalous super-diffusive regime found in active 
bacterial suspensions . Other regimes include inertial 

or propulsive behaviors, diffusive behavior ~ and 
an anomalous correction of the form —jt^^^. 

We consider a spherical colloidal particle of radius R 
with an axially symmetric pattern of surface activity as 
shown in Fig. [^h-, which in essence leads to the re- 
lease of (excess) product particles V with diffusion co- 
efficient D. This could correspond to a chemical reaction 
S — > S' +V catalyzed on the surface of the colloid, with 
the simplifying assumption that the produced iS' will act 
almost like the consumed S (see Fig. O)). On the other 
hand, it could also correspond to a system that actively 
releases P particles from the interior of the sphere. The 
product particles will have a concentration profile C(r, t) 
that could interact with the fluid in the vicinity of the 
colloid due to surface phoretic effects, and cause relative 
motion with a slip velocity 'Vs{9,(j),t) = iJ,\7\\C{R,9,(t),t) 
(in standard spherical coordinates; sec Fig. [2^), where fi 
is the surface mobility and Vy denote the lateral gra- 
dient In diffusiophoresis, fj, = kBTX^/r], where 
is the thermal energy scale, 77 is the viscosity of 
water, and A is the Derjaguin length [iJl, which is de- 
fined in terms of the interaction potential W(z) between 
the diffusing particles and the surface of the colloid as 
A2 = dz z[l- e-W'(^)/fcBTj _ j^Q^g ^i^g^^ Q^j. 

tions, will be positive (negative) for repulsive (attrac- 
tive) effective surface potentials. In the case of ionic sys- 
tems A will be set by the Debye length d, [l^l . Averaging 
over the surface of the sphere (using integration over the 
solid angle Q), we can obtain the instantaneous velocity 
of the colloid as v(t) ^ J dn /iV||C(i?, 6*, </>, t). 

The axis of symmetry of the colloid, which points 
to the direction of propulsion for sufficiently asymmet- 
ric patterns 0], is defined by the unit vector n{t) = 
(sin 9n (t) cos (fin (t) , sin 0„ {t) sin </>„ (t ) , cos 9n{t)) (sec [5^) . 
The stochastic nature of n(i) due to rotational diffusion 



of the colloid causes the cloud of product particles to con- 
stantly redistribute, which will in turn make the velocity 
of the active colloid fluctuate. To get the instantaneous 
velocity of the colloid, we need to solve the diffusion equa- 
tion for the concentration profile of the product particles, 
namely 

dtC{v, t) - DV^C'ir, t) = a{9, 0, t)5{r - R), (1) 

subject to the boundary condition of vanishing normal 
current on the surface of the sphere. In Eq. ([T]), a{9, 0, t) 
is the surface activity function of the sphere, i.e. rate per 
unit area of the introduction of (excess) product parti- 
cles. For axially symmetric surface activity, we can repre- 
sent the function in terms of the spherical harmonics as 

a{9,<i>,t) = Y.L.^{^)<^t Y;Mt),Mt))Yim{0,^). 

Equation ^ only gives the average density, and the lin- 
ear relation between the velocity and the concentration 
profile suggests that in order to calculate velocity corre- 
lations we need to incorporate the density fluctuations as 
well, which we do using the method outlined in Ref. 0. 

Using the formulation described above we can cal- 
culate the velocity autocorrelation function for the ac- 
tive colloid A"^{t) = (v(t) • v(0)) as well as the mean- 
squared displacement AL^{t) = ([r(t) — r(0)]^) = 

dti d<2 (v(<i) • v(t2)) . There are three important 
time scales in the problem. The characteristic diffu- 
sion time of the product particles around the sphere is 
Td = R^/D, where D = ksT / {Gnrja) depends on the ra- 
dius of the solute particles a. This time scale sets the re- 
laxation time of the redistribution of the particles around 
the sphere when it changes orientation. The rotational 
diffusion time, Tr = ^TTrjR^ / k^T ^ controls the changes in 
the orientation of the sphere, and is defined via the ori- 
entation autocorrelation function: (n(t) • n(0)) = e"*/'^'' 
[ist . Finally, the hydrodynamic time that controls the 
crossover between the inertial and the viscous regimes 
is given as Th = R? jv^ where v = rj/ p is the kinematic 
viscosity of water that involves the mass density p. Prac- 
tically speaking, we always have <^ Tr although 
this is not a fundamental requirement. 

We can identify three distinct contributions to the ve- 
locity autocorrelation function: (1) a contribution from 
the density fluctuations, which turns out to be only sen- 
sitive to the overall symmetric component of the activity 
and is present even for non-propelled active colloids, (2) 
a contribution from the asymmetric component of the ac- 
tivity, and (3) a hydrodynamic contribution that entails 
the passive diffusion of the colloid and the hydrodynamic 
long-time tail. These will lead to distinct contributions 
to the mean-squared displacement, which will add up 
to make the total mean-squared displacement, namely 
AL^t) = AL%^{t) + ALLytn W + ^L^^t). We wiU fo- 
cus on each of these contributions separately below. 

Symmetric Contribution. — The density fluctuations 
that are accounted for by adding a noise term to Eq. 
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([T|) lead to a contribution to the velocity autocorrelation 
function that is proportional to ao [ll] — the i = coeffi- 
cient in the expression for a{9, </), t) in terms of the spher- 
ical harmonics [ag = / dft a{9, </>, t)]. This means that 
the contribution by density fluctuations is only sensitive 
to the mean overall surface activity of the colloid, and 
not the patterning structures on it. We find the asymp- 

;^^§^(VTrf)-i/2 for t « r, 

32|aw(V^<i)"'/' for t » Td. Wc also 
calculate the contribution of density fiuctuations to the 
mean-squared displacement, which has the asymptotic 
behaviors of 

t'/' ;t«r,, (2) 
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at long times, and a smooth crossover between them. 
Here, ci = 1.17810 is a numerical prcfactor. 

Asymmetric Contribution. — Solving Eq. ((T|) with- 
out the noise term, we can calculate the propulsion ve- 
locity of the colloid as a function of time for a given 
time dependent orientation trajectory. We find v(t) = 
^ dt'M{t - t')n{t'), where vq = ~ai^/{3D) is the 
mean propulsion velocity and the memory kernel is 
given as M{t) = ^ /o°° c?"^^T^e~"'^*/^'^\ with asymp- 
totic behaviors M(t) ~ -^[t/rd)^^^'^ for t 'ti and 
M(t) ~ g^(i/rd)-5/2 for t > Td. Note that the propul- 
sion velocity is controlled by the £ = 1 term (ai) in the 
surface activity profile. 

Rotational diffusion of the colloid randomizes its ori- 
entation over the time scale r^, which leads to a con- 
tribution to the velocity autocorrelation function of the 
form of a convolution between two memory kernels and 
the orientation autocorrelation function. This leads to a 
velocity autocorrelation function [l6| . which have three 
different regimes, due to the presence of two characteris- 
tic time scales Td and Tr- At short times, t Td Tr, wc 
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where C2 = 0.642699 is a numerical prcfactor, which im- 
plies that the autocorrelation function will be rounded 
off at small t. For intermediate time, <C t ^ 
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for long times, Td <C Tr 

2 L-t/r, I _3_ 3/2,_5/2 
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cay at long times is primarily algebraic and not expo- 
nential. Consequently, the mean-squared displacement 
will have three different regimes. We find the asymptotic 
form of 
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at short times. 
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at long times, with a smooth crossover between them. 

Hydrodynamic Contribution. — Thermal fluctuations 
of the solvent fluid velocity also contribute to the velocity 
autocorrelation function of the sphere, because of the no- 
slip boundary condition between the fluid and the colloid. 
Performing a similar calculation, one finds [l3| AL^^^j ~ 

t'/' for t » Th, where Do = kj^T / {(Sttt^R) 
is the bare diffusion coefficient of the colloid. The first 
term in the above equation describes the standard passive 
diffusion of the sphere while the second term corresponds 
to the hydrodynamic long-time tail [3, US]- At short 
times when t <C r/i, one finds AL^^^ ~ 3 (y^MTt) where 
Moff is the effective inertial mass of the colloid in water. 
The above results are summarized in Fig. [T] 

Discussion. — At the longest time scales [t > Tr), 
all of the contributions arc diffusive, leading to a total 
effective diffusion coefficient 



kBT inaffj.'^'nR^ 



(7) 



The different terms in the above expression exhibit dif- 
ferent /^-dependencies, which causes the asymmetric con- 

1/2 

tribution to be dominant for R > [DksT / (aiiirj)] , 
while the symmetric contribution takes over when R < 
ctQ^i^rj^ /{D'^k-eT). At the shortest time scales, on the 
other hand, the contribution due to phoretic effects will 
also be dominated by inertial effects that should lead to 
baUistic contributions (see Fig. [T|). 

In the intermediate times, we observe a number of 
anomalous behaviors. For symmetric active colloids, the 
super-diffusive t^^"^ behavior [Eq. ([2])] for Th < t < Td is 
a new regime (for isolated self-propelled particles) where 
the motion is neither ballistic nor diffusive. The rea- 
son a symmetric particle can move at all is because den- 
sity fluctuations of the cloud of solute particles can in- 
stantaneously produce an asymmetric distribution and 
therefore a net propulsion in some direction. This mo- 
tion, however, will be decorrelated via density fluctua- 
tions themselves, leading to fiuctuations without sym- 
metry breaking. We can understand the form of Eq. 
([2]) as follows: using AL^ ~ v'^t'^, and putting v ^ 
fiVC ~ ixSC/R, we find AL^ - fjL'^{SC{t)SC{0))t'^ /R^. 
The density auto-correlation function can be written as 
{5C{t)dCiO)) = (SC^)k{t), involving the density fluctua- 
tions ((5C^) and the kernel k{t) that controls the relevant 
relaxation mode (20| . Here, relaxation is controlled by 
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diffusion, hence k(t) ^ l/{Dt)^/^, and the number fluctu- 
ations are controlled by the average number of particles 
{{SN"^) ^ -/Vave) — as inherent to any Poisson process — 
that yield {SC"^) ~ Cave- On the other hand, the average 
density is controlled by the average particle production 
rate (per unit area) ao as Cave ^ {aoB?t) / . Putting 
these all together, we find Eq. ([2]). This shows that the 
active velocity fluctuations are controlled by two mech- 
anisms: particle production (that controls the density 
fluctuations) and diffusion of the produced particles. In- 
terestingly, a similar t^^-^ power law has been observed in 
the motion of passive tracer particles in a bath of bac- 
teria whose flagella stir up the fluid (and theoretically 
accounted for using aphenomenological continuum ac- 
tive medium theory) [13] , and in ion channel gating [2l[ , 
both of which cases are also governed by some sort of 
density fluctuations p^ . 

For asymmetric particles when < t < Td, the t'^/^ 
contribution is added (with a positive coefficient) to the 

propulsive term. On the other hand, for Td < t < 
the memory effect that exists for self-propelled asym- 
metric colloids introduces an anomalous anti- correlation 
(i.e. contribution with negative sign) in the velocity auto- 
correlation function and the mean-squared displacement 
[Eq. ([5])]. Such anomalous corrections, which have also 
been observed in continuum theories of interacting active 
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self-propelled particles [ll| , are reminiscent of the effect 
of the hydrodynamic long-time tail. Note, however, that 
the anomalous —"ft'^/'^ correction in Eq. ([5]) corresponds 
to much longer time scales and should be more easily 
observable than the hydrodynamic long-time tail. 

To get a better feel for the working domain of each 
regime, we can write the time scales (for water at room 
temperature and using a typical value of a = 1 A) in 
the following convenient form: = 10^^ (-R/l/im)^ s, 
Td = 10^^ {R/liiu\f s, and = 3 {R/liim)^ s. This 
shows that while rapid cameras or scattering techniques 



19j could in principle resolve all the three domains 
for micron-sized beads, using i? = 20 /im (which yields 
Td = 0.4 s and = 2.4 x 10^ s) should provide a comfort- 
able working range for an experiment that aims to resolve 
the anomalous components of the motion. Finally, we 
note that here wc have only focused on the phoretic con- 
tributions to velocity fluctuations, and in practice other 
sources of fluctuations might also be present [lo| , which 
need to be taken into account. 

In summary, we have shown that active colloidal parti- 
cles that interact with their self-generated cloud of solute 
particles can have a range of different types of stochastic 
motions. Using parameters such as surface activity, sur- 
face mobility, and size, we can tune the behavior of active 
colloidal particles, and this could provide new possibili- 
ties in designing functional motile agents for applications 
in micro- and nano-fluidics and targeted delivery. 

I acknowledge invaluable discussions with A. Ajdari, 
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